Asymmetric binaries and hierarchy of states in self-gravitating systems 
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We complete the analysis of the equilibrium shapes of rotating self-gravitating gases initiated in 
Phys. Rev. Lett. 89 031101 (2002) by studying the formation and thermodynamics of asymmetric 
binaries and by introducing a family of order parameters to classify the possible equilibrium states. 
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Gravity-dominated systems such as classical self- 
gravitating gases are known to possess highly non- 
standard thermodynamic properties ||, 1). The long- 
range nature of the Newtonian potential prevents ther- 
modynamic functionals from being extensive Q|, caus- 
ing such features as the inequivalence of statistical en- 
sembles, the occurrence of negative specific heats and 
of collapse transitions, and the existence inhomogeneous 
ground states ||, ^, |^, ^, These attributes are not 
peculiar to gravitating gases alone. Similar traits mark 
the behaviour of systems with long-range, non-integrable 
potentials [0, of finite systems or of systems of a size 
comparable to the interaction range, such as atomic clus- 
ters or highly excited ("hot") nuclei 
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are even 



present in certain short-ranged models [13| and in general 
at phase separations [Q . Besides the relevance to astro- 
physics and cosmology strictly speaking, understanding 
the equilibrium state of self-gravitating systems is hence 
important also from a broader, fundamental viewpoint 

In a recent work |l^, |l^, to be later referred to as I, 
we investigated the possible equilibrium shapes of three- 
dimensional rotating self-gravitating gases using a mi- 
crocanonical mean-field theory in which the system is 
enclosed in a finite spherical volume V with the total en- 
ergy E and the total angular momentum L conserved. 
The singularity of the Newtonian potential at short dis- 
tances was cured by the use of Lynden-Bell statistics [|| . 
It turned out that, along with the known "gaseous" and 
"collapsed" (a single, dense cluster) states, the system 
can form symmetric "binaries" (two identical clusters) if 
E is sufficiently low, that is if gravity is strong enough, 
and rotation is sufficiently fast. Moreover, at intermedi- 
ate E and L, a large phase-coexistence region with neg- 
ative specific heat appears, where the system is found 
in dense structures (either a single cluster or a binary) 
embedded in a vapor. 

The purpose of this paper is to complete the study 
of I by describing the formation and thermodynamics of 
asymmetric binaries (which are far more common in the 
universe than symmetric ones). Furthermore, a family of 
order parameters will be introduced to discern the pos- 
sible equilibria and allow for a classification of the equi- 
librium states. Some "exotic" equilibria (with a very low 
entropy) that can be obtained using these order param- 
eters are presented as an example. 

We begin by recalling the key steps in the theory of I. 



The Hamiltonian is = 1, 



H = 



(1) 



with Ti^VdW" denoting the position and S 
the momentum of the i-th particle, to being its mass. 
G is the gravitational constant. The task is to find the 
particle density profiles that maximize the Boltzmann 
entropy S — k log W, where the state sum W is given by 



W : 
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J 5{H - E) 5{L - ^ X d^^p (2) 



with a a constant. is a function of E and L, as are 
S and the desired density profiles. The radius R of the 
box is taken to be fixed. Adhering to the notation of I, 
we shift to the reduced variable x = r/R, which we shall 
indicate as a; = (xi, X2, x^) (Cartesian representation) 
or X = {x,6,(f>) (spherical representation). Particles are 
assumed to have "hard cores" , so that the density profile 
c is normalized as / c{x)dx — 8, with 0<c<l. 0isa 
fixed parameter (0 < < 47r/3) which throughout this 
study, as in I, is taken to be = 0.02. A brief discussion 
about the dependence of the results on Q is given at the 
end of the paper. Finally, we take to = 1 and measure 
energy and angular momentum in units of GN^/R and 
V RGN^, respectively. 

Assuming Lynden-Bell statistics to mimic the hard 
cores Q , at stationary points of the entropy surface the 
density profile c satisfies the condition (see I for details) 
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/3 = 



^^ dx' + ^/3{ujxxr-t, (3) 
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[E - \L^l-^L - $[c]] ' 



(4) 



u) = I~^L is the angular velocity, and /x is a Lagrange 
multiplier implementing the constraint on 0. I stands 
for the inertia tensor, with elements (a, 5=1,2, 3) 
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c{x){x'^6ab - XaXb)dx 



(5) 



(in units oi NR^), while <I> is the gravitational potential. 
In general, the solutions of Eq. (0), whose number and 
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form depend on E and Zi, can be written in spherical 
coordinates as 



(6) 



1=0 m=-l 



where the "weights" bim{x) have to be determined by 
solving (numerically) the system of integral equations 



= / g(x,e,(j))Yun{e,<i>) dcos0 d0 (7) 



1 + !:,,,„ u,™(:i:)yi™(0,0)-^/3c^'a;"sin2e+^; 



Some details of this task are discussed in I. 

If L is taken to lie along the 3-axis and if only even 
harmonics (/ — 0, 2, 4, . . .) are included in the series (|^), 
entropy-maximizing solutions of (|^) can be divided in 
two main classes (see I): (a) axial-rotationally symmetric 
ones, e.g. gas-cloud (high i?, low L), distorted gas clouds 
or disks (high E, high L), and single clusters (low E, low 
L); (b) axial-rotationally asymmetric ones, like identical 
double clusters or symmetric binaries (SBs, low E, high 
L). The former are independent of the 4> angle, and thus 
for them the series (^ involves only the zonal harmonics, 
with rn = 0. The latter, instead, must depend on 
and the lowest-order term in (^) that can provide this 
dependence is the sectoral harmonics i2,2- In fact, on 
physical grounds, SBs can be characterized by a non-zero 
value of the order parameter 



(8) 



which is proportional to the difference between the di- 
agonal components of the inertia tensor in the 1 and 2 
directions. Using (||) and (||) it is indeed easy to see that 
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(a) B: 



-0.60, L = 0.46 



(b) E = -1.16, L = 0.46 



FIG. 1: Two asymmetric binaries obtained from Eq. (|3|). 
The contour plot is shown together with the density profile, 
(a) Incipient binary, with a broad low-density body detached 



from the central cluster, 
energy). 



(b) Well-formed binary (very low 



energy. As a consequence, for SBs the "random" kinetic 
contribution E — L^/(2/) is larger than for ABs. Thus 
they have a larger entropy. So ABs are expected to be 
strongly suppressed or significantly less probable at high 
angular momenta. Indeed, we were unable to find stable 
asymmetric binaries for L > 0.5. The fact that ABs tend 
to form the smaller cluster in proximity of the bound- 
ary of the box is due to the fact that ABs start to form 
by evaporating small masses from a large, dense cluster. 
The conservation of the center of mass relegates them to 
the border. In well formed binaries, the peak of density 
in the smaller cluster is separated from the boundary. 

In Fig. H the caloric curve f3 = dS/dE versus E is 
shown at L = 0.46. The jump occurring at E = —0.2, 
where ABs become more favorable than SBs, is an arti- 
fact of the asymptotics employed in the mean-field theory 
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(9) 



The inclusion of odd harmonics leads to new solutions, 
namely asymmetric binaries (ABs). A sample of asym- 
metric binaries obtained in this way (with the addition 
of fixing the center of mass at a; = 0) is shown in Fig. ||. 
As it was to be expected, such solutions appear at low 
energies at intermediate angular momenta (L ~ 0.46) 
and exist in a small range of values of L up to L ~ 0.5. 
This whole region lies in the mixed phase with negative 
specific heat of I, hence asymmetric binaries at equilib- 
rium are embedded in a vapor and they don't exist as a 
"pure" phase. Moreover, in this range they are signifi- 
cantly more probable than symmetric binaries, as shown 
in Fig. H where the entropy of the different stable solu- 
tions at L = 0.46 is shown. 

A simple physical argument suffices to explain why at 
high angular momenta ABs must have a much lower en- 
tropy than SBs, if they exist at all. SBs have a larger 
moment of inertia than ABs, hence a smaller rotational 
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FIG. 2: Entropy of the different stable solutions of Eq. (y) 
at L = 0.46. The quantity has been subtracted from 
the entropy to make the separation and the curvature more 
evident. 
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FIG. 3: Caloric curve (3 vs E a.t L = 0.46. Notice that the 
curvature of the entropy surface as a function of E and L, i.e. 
SeeSll ~ S%i^, is negative along the AB branch. 



of I. (Being S a smooth and multiply-differentiable func- 
tion, such a feature is clearly spurious.) 

Again on physical grounds, the natural order parame- 
ter for AB solutions is (xf ) := / x\ c{x)dx, because in 
spatially asymmetric structures it must be non-zero be- 
cause of evident symmetry reasons (the two clusters are 
aligned along the xi axis), at odds with what must hap- 
pen in symmetric structures. Now using (^) one can see 
by working out explicitly the integrals involving spherical 
harmonics |17|| that 
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x^b3^i{x)dx- 



x^b3^3{x)dx (10) 



In principle, all of the above terms should be considered. 
However, in any spatially symmetric structure (gas, sin- 
gle cluster, SBs) all of them are zero, because odd har- 
monics do not contribute, as discussed above. Each of 
them is non-zero only in presence of spatial asymmetries, 
as in ABs, hence any of them can be used as an order 
parameter. In particular, we find it convenient to choose 
the quantity 



-03,1 = 



/ 



x^ b3^i{x)dx 



(11) 



whose behaviour at L = 0.46 is shown in Fig. |j. AB 
solutions appear at E = —0.18, where the order parame- 
ter jumps discontinuously from zero to a finite value, and 
becomes more favorable than SBs aX E — —0.2. 

Eq.s (H) and ( [III ) suggest that, in general, one can 
classify the different equilibrium states by introducing 
the family of order parameters defined by 



Di^m — 



x'-^^bi^rn{x)dx 



(12) 



with integers I = 0, 1, 2, . . . and m = —I, + 

For small I, each of these functions has a simple physical 

interpretation. 

For I — 0, -Do,o is simply a constant proportional to 0, 
since, recalling that / c{x)dx = (1) = Q, 



z(x)dx — j x'^bQ^o{x)d2 



(13) 



For / = 1, one has that -Di.i, Di q and i^i,-! are noth- 
ing but the coordinates of the center of mass (modulo a 
proportionality constant) : 



„CM 



■^1 




x^bi^i{x)dx 



x^'^' = {X2) = \IyJ ^'K~i(.^)dx (14) 

X^'' EE (X3> - / ^'h,0ix)dx 

For I — 2, the situation is slightly more involved. 152,0 
can be easily seen to be linked to the quantity 



{x' - 3x1) 
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x^b2fi{x)dx 



(15) 



which is a convenient order parameter for disk (deformed 
gas) solutions. In fact, for an hypothetical configuration 
being a perfect sphere (say a homogeneous gas cloud fill- 
ing the whole V), one would have I?2,o = evidently. 
When rotation deforms the cloud by shrinking it along 
the 3-axis, one would have D2.0 7^ 0. As discussed in 
I, is instead the order parameter discerning axial- 
rotationally symmetric structures from axial-rotationally 
asymmetric ones (it is non zero for both SBs and ABs). 
It is particularly convenient for SBs because of its con- 
nection to the inertia tensor components, see (|9|). As for 
the remaining -D2,m's, they can be easily seen to be re- 
lated to the off-diagonal components of the inertia tensor 
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FIG. 4: Behaviour of the order parameter for asymmetric 
binaries 1)3,1 as a function of energy at L = 0.46. 
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FIG. 5: Two "exotic" solutions of (g) obtained via the order 
parameters (p^). (a) Single cluster surrounded by a ring, (b) 
Single cluster with two "satellites". These solutions coexist 
at E ^ -0.5 and L = 0.7. 

It is clear that higher-order D's are all connected to 
homogeneous polynomials in xi, X2, x^, which are in turn 
connected to the Cartesian representation of spherical 
harmonics For example, for Z = 3, D^ q is connected 
to 



(2x1 - 3xjx3 - 3x1x3) = ]/— J x^b3,oix)dx (16) 

Their physical interpretation for larger I is however less 
straightforward. 

Using these order parameters, solutions can be classi- 
fied hierarchically, a procedure that is particularly effec- 
tive at high L. In our setting, with the center of mass 
fixed at a; = one has the following picture. Homo- 
geneous gas clouds have all Dim's equal to zero except 
Do Q. For deformed clouds (disks) only Di q with I even 
are non zero. For symmetric binaries (which are con- 
nected to "gas" -like solutions continuously, see I), also 
^2,2 7^ (and in general all Di^„i with / and m even), all 
other D's being zero, while for asymmetric binaries -Da,! 



is also not zero. This classification can be continued fur- 
ther, for the order parameters (jj) allow to construct all 
possible solutions of (||), by adding a proper field that 
selects the solution for which the chosen order parame- 
ter is not zero. As an example, we show in Fig. || two 
of several equilibria that can be obtained in this way at 
very high angular momentum. Of course, such "exotic" 
solutions of (|3|), if stable, have much lower entropy than 
those discussed in I and here. 

To summarize, we have carried further the equilibrium 
study of rotating self-gravitating gases of classical non- 
overlapping particles started in I by describing the forma- 
tion and thermodynamics of asymmetric binaries. More- 
over, we introduced a family of order parameters that 
allow to characterize all the different shapes that can 
be encountered (homogeneous gas, disk, symmetric bi- 
nary, asymmetric binary, along with more exotic states) , 
and discussed their physical meaning in simple cases. 
The present study does not substantially alter the global 
phase diagram presented in I, since asymmetric binaries 
do not exist as a pure thermodynamic state but only in 
the phase coexistence region. Among the problems that 
remain open we would like to mention the following two. 
First, the dependence on 6. Our results have been ob- 
tained at fixed Q = 0.02. Different O's may lead to dif- 
ferent results. By letting 8 — > one obtains the case 
of overlapping particles (hard cores are neglected and 
Lynden-Bell statistics is replaced by Boltzmann statis- 
tics in space), which is a classical problem in cosmology 
dating back to The rotating case is dealt with in 

|l9| . For Boltzmann particles we haven't found any evi- 
dence of asymmetric binaries. On the other hand, when 
Q grows (O < 47r/3) the system may become too closely 
packed and the interesting features one observes for small 
O might disappear at some point. The second problem 
is the dependence on R. Throughout this paper we have 
assumed that R is fixed but it would be interesting to 
analyze this problem in a context in which R varies, for 
instance to mimic an expanding-universe scenario. 
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